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Abstract 


The Comment by Aifantis that criticizes the article ‘On non-singular crack fields in Helmholtz type enriched elasticity 
1 —' theories’ [Lazar, M., Polyzos, D., 2014. Int. J. Solids Struct, doi: 10.1016/j.ijsolstr.2014.01.002] is refuted by means of 
^ clear and straightforward arguments. Important theoretical aspects of gradient enriched elasticity theories which emerge 
I in this work are also discussed. 
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A comment by lAifantisI ( 20141) [thereafter referred 
to as (A)l criticizes som e aspects of the arti¬ 
cle ( Lazar and Polvzosl 2014 1 [referred to as (LP)]. In 
the (LP) article, we discuss and point out the physical 
or unphysical meaning and interpretation of simple non¬ 
singular s tress fiel d s of c r acks o f mode I and mode III pub- 
lis hed by Aifantis 120091 201ll); Isaksson et al.l 12012 1 and 
llsaksson and Hagglundl (12013 1. In particular, we investi¬ 
gate if the aforementioned solutions satisfy the equilibrium 
conditions, boundary conditions and compatibility condi¬ 
tions in the framework of nonloca l or gradient elasticity 
theory. I n the meantime, Konstantonoulos and Aifantid 


tneory . i n tne meantime. iKon stantopoulos and Ailantis 
(2013) and Isaksson and Dumontl (l2014h have continued to 


X 


publish their non-singular stress fields in other articles in 
scientific journals. T he aim o f the (LP) article was to show 
that the results of (lAifantid. 

120121 : Isaksson and Hagglundl 


2009. 2011; Isaksson et al 


20131) cannot be the correct 


solutions of a nonlocal or gradient compatible elastic frac¬ 
ture mechanics problem, since important physical condi¬ 
tions are violated and not to discuss about the physical 
assumptions of the GRADELA model made by Aifantis 
and others. In this work, we reply to the comments of 
(A) which are relevant to the (LP) article and we clarify 
important aspects of gradient enriched elasticity theories. 
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2. Discussion 


We review in the (LP) ar ticle the f r amework of Eringen’s 
nonlocal elasticity theory dEringeil ll983[ 12002 ). In this 
context, (A) claimed that “Eringen never uses the classical 
stress cr b in a balance lav0 a ifj = 0” > writing additionally 
that “Eq. (LP5) is an artifact of a number of assumptions 
and not a fundamental equation of Eringen’s theory”. It 
is easy to show that (A) is mistaken concerning this point 
due to the following reasons. First of all, if the tensor cr° 
is the classical stress tensor, then it has to fulfill certainly 
the classical equilibrium condition oA ■ = 0. Moreover, 
combining Eq. (LP4) with Eq. (LP7), one finds 

= a ij,j = 0 ) (1) 

which proves that the relation cA ■ = 0 is true and is part 
of Eringen’s theory of nonlocal elasticity. Here, 

L=l-l 2 A (2) 


is the characteristic differential operator appearing in Erin¬ 
gen’s nonlocal elasticity and in the corresponding gradient 
elasticity, tij is the stress tensor of nonlocal elasticity, £ 
is a characteristic length scale and A denotes the Lapla- 
cian. L is usually called as Helmhol tz operator in gra- 


(1995 

); 

); 

Vardoulakis et al. (19961): Gutkin and Aifantis 

(1999 

Peerlings et al. (2001): Paulino et al. ( 

2003); 

Gutkin and Ovid’ko (20041): Lazar and Man gin ( 

200|); 


1 Actually, oA . = 0 is a conservation law. 
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Chan et al. ( 200fih : Aifantisl ( 20091) : Lazar ( 2013 )). In gen¬ 
eral, the Helmholtz equation Am + cm = 0, c = constant, 
can be found with positive (c > 0) or negative sign (c < 


0) in th e mat h emat ic al literatu r e (e.g ., Koshlvakov et all 


( 19641) : IroosI (Il969l) : IPolvaninl rt200ll) : iHassaxiil ri2002l) ): 


and only sometimes the Helmholtz equation with nega¬ 
tive sign, that is when c < 0, is referred to as the mod¬ 
ified Helmhol t z equation (e.g., IZauderen (|1983l) ). More¬ 
over, Eringen ( 19831 2002 ) stated that in the static case of 
nonlocal elasticity, the classical equilibrium condi tion (Eq. 
(LP5)) is to be satisfied (see, e.g., Eq. (6.9.32) in Eringenl 
(20Q2j)). Thus, the classical equilibrium condition is part 


of Eringen’s theory of nonlocal elasticity in contrast to the 
claim of (A). In addition, we note that anisotropic nonlo¬ 
ca l elasticity of generaliz ed Helmholtz type can be found 
in Lazar and Agiasofitou (1 2011 ). 


Afterwards, we review in (LP) a simplified gradi¬ 
ent elasticity theory, the so-called gradient elasticity 
of Helmholtz type. This simplified and straightfor¬ 
ward version of gradient elasticity with only one length 
scale i is a particular case of Mi ndlin ’s th eory of gra¬ 
dient elasticity of form II ( Mindlinl . Il964 ) and it has 
been called bv lLazar et al. ( 20051) “gradient elasticity of 
Helmholtz type”. The development of grad ient elasticity 
of Helmholtz type ( Lazar and Mauginl . 120051 ) was strongly 
influenced and stimulated by the field theoretical frame¬ 
work of the so -called dislocation gauge theory ( Lazail 
2002a b. 2003allt]) and by Feyn man’s lectures on gravita- 
tion (Feynman. 19951) (see also Lazar and Anastassiadisl 


(2009 ): lAgiasofi tou and La,zan (120101)). As shown 


Lazar and MauginT i 2005 ): Lazar ( 2013 . 20141) . gradi¬ 


ent elasticity of Helmholtz type is based on proper vari¬ 
ational principles including double stresses and serves 
a “physical” regularization ba sed on higher-order par¬ 
tial d if ferenti al equa t ions. In Lazar and Mauginl ( 20051 . 
20061) : L aza r J20ld 120131 l2014h gradient elasticity of 
Helmholtz type has been successfully applied to calcu¬ 
late non-singular stresses, non-singular strains and non¬ 
singular displacement fields of straig ht dis lo cations and 
disloc ation loops. It is noted that lLazar and Maugiri 
( 20051) is the first work where such a task was under¬ 
taken within the framework of Mindlin’s gradient elas¬ 
ticity with the remarkable outcome that both the stress 
and the strain singularities are removed at the dislocation 
line in contrast to th e corresponding solutions of classical 
elasticity. Moreover, lLazar and Maugin ( 20061) corrected 
the mistaken gradient elasticity solutions for the displace- 
me nt fields of strai ght sc rew and edge disl ocations given 
by Gutkin and Aifantisl ( 19961 ll997 . 19991) . Later, (A) 
adopted the correct “Lazar-M a ugin solutions” in further 
studie s ( Davoudi and Aifantisl . 12013 ) (see, e.g., Gutkinl 
( 20061) : Davoudi et al.1 ( 20ld )). It should be mentioned 
that the non-singular dislocation key-formu las obtained 
in gradient elasticity of Helmholtz type (lLazarl . 2013, 
20141) have been rece ntly impl e ment ed in 3D discrete 


dislocation dy namics ( Po et al., 2014[). The J-integral 


ble gradient elasticity of M indlin and Helmholtz type and 
by lAgiasofitou and Lazar ( 20091) for Mi ndlin ’s g radient 
elasticity of grade three. In addition. lLazar ( 2014 ) shows 
that the framework of gradient elasticity of Helmholtz type 
is in agreement with the framework of the magnetostatic 
gradient theory, which is the magnetostatic part of the so- 
called Bopp-Podolsky theory ( Bum ). 1940 ) Podolsky , 194 20 
and how gradient theories are used in_physics. However, 
it has to be mentioned that iMindlin ( 1964 ) introduced his 
theory of gradient e lastici ty without giving any credit to 
Bopp and Podolsky. Gunther ( 1976 ) was the first to speak 
of a mechanical model of the Bopp-Podolsky potential for 
defects in elasticity. 


On the other h an d, abo ut three decades after Mindlin 
(Il964l). Aifantisl (1 19921): Altan and Aifantisl ( 1992 ) 


R.u and Aifantisl dl993l) : Gutkin and Aifantid ( 1996 . 1997 ) 
introduced a simp le gradie nt elasticity model called 
GRADELA. In fact, Aifantisl ( 19921 ) has intuitively postu¬ 
lated a gradie nt modifi c ation of Hooke’s law. However, th e 
framewor k of Aifanth ( 1992 ): Altan and Aifantisl ( 1992 ): 
Gutkin and Aifantisl ( 1996 Il997l) lacks double stresses 
and, consequently, the stress tensor remains singular. 
Their approach postulates a sophisticated Hooke’s law 
including a Laplacian of the elasti c strain tensor. In 

Gutkin and Aifantis 


a refined version of their model. 


( 19991) found non-singular stress and non-singular strain 
fields using a more sophisticated Hooke’s law including 
a Laplacian of the (Cauchy) stress tensor. Moreover, 
in c ontrast t o Mi ndlin’s the ory of gra d ient elastic- 
ity ( Mindlin . 19641) ( s ee al s o Ijaunzemls ( 1967J) ) , the 
fra mework of Aifantisl (f l992h: Al tan and Aifantisl (1992) 


and Gutkin and Aifantisl ( 1996 . 19971 1999fT is not based 
on proper variational considerations (e.g. to obtain per¬ 
tinent boundary conditions). Thus, GRADELA is based 
on ad-hoc postulations and not on a prop er theoretical 


frame work as basis of a theory (see also [Fafalis et al 


( 2012 ) for an interesting discussion). It becomes evident 
that gradient elasticity of Helmholtz type is different 
from GRADELA due to the aforementioned reasons and 
that the authors of the (LP) article have not adopted 
GRADELA. 

Now, we sketch t he basics of Min dlin’s theory of gradient 
elasticity of form II ( Mindlinl . 11964 ) and we give its relation 
to the theory of gradient elasticity of Helmholtz type. In 
Mindlin’s theory of gradient elasticity of form II, the strain 
energy density has the general form 


— W (&ij , ‘ 


(3) 


where eij denotes the elastic strain tensor of gradient elas¬ 
ticity. T he stresses of gradie nt elasticity are defined as 
(see, also Eshel and Rosenfeldl ( 197(1 1973 )) 


dW 


Tij ~ de T ' 


ji 


was given by lLazar and Kirchner ( 2007 ) for incompati 


dW 

1~ijk — — 1~jik 


ij,k 


Cauchy stress tensor , (4) 

double stress tensor. (5) 
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Thus, Tij is the Cauchy stress tensor of gradient elasticity 
which usually reads 


~/j — CijklCkl 


( 6 ) 


where Cijki is the tensor of the elastic constants. In the 
(LP) article, we use the notation Cauchy-likc stress tensor 
for , which means Cauchy stress tensor of gradient elas¬ 
ticity. The eq uilibrium condition for the stresses takes the 
form ( Mindlinl . 1964h 


^~ijk : kj — 0 • 


(7) 


The relation of Mindlin’s theory of gradient elasticity of 
form II to gradient elasticity of Helmholtz type is given by 
the particular form of the double stress tensor © 

Tijk — ^ Tij^k • ( 8 ) 


The double stress tensor © is nothing but the gradient of 
the Cauchy stress tensor of gradient elasticity multiplied 
by the pre-fact or l 2 in order to have the c orrect dimen¬ 
sion ( see, e.g., Lazar and Mauginl ( 2005h : IPolvzos et al. 


(20031)). Eq. © can be easily obtained by Eq. © us¬ 
ing the strain energy density o f Mind lin’s theory of Form 
II (see Eq. (11.3) in Mindlinl ( 19641 ) 1 with an appropri¬ 
ate simplification of the five gradient coefficients (see, 
e.g., lLazar and Maugiiil ( 20051) 1. The equilibrium condi¬ 
tion © with the help of the relation © can be reduced 
to 


-T' T ij,j — 0 , (9) 

where L is the Helmholtz operator given in Eq. ©. 
Eq. © is the equilibrium condition in gradient elastic¬ 
ity of Helmholtz type. Furthermore, substituting Eq. © 
and a compatible elastic strain tensor eij = 1/2 (uij + Uj^) 
into Eq. ©, one can find a homogenous Helmholtz-Navier 
equation for the displacement vector Uk 


LCijkiUk,ij = 0 , ( 10 ) 


which is an elliptic partial differential equation of fourth- 
order. 

Using an operator spl i t (see , also Vekual ( 1967 1: Lazar 
( 20141) : Ir.u and Aifantid ll99,^ 1. one can decompose the 
partial differential equation of third-order © into a sys¬ 
tem of partial differential equations consisting of a partial 
differential equation of first-order and a partial differential 
equation of second-order, as follows 


a ij,j = 0 . 
LTij = . 


( 11 ) 

( 12 ) 


Eq. © is an inhomogeneous Helmholtz equation where 
the inhomogeneous part is given by <r®., which has to fulfill 
Eq. ED. Since cr/ ; satisfies the classical equilibrium con¬ 
dition (fill) and does not depend on the length scale £, it 
may be identified with the Cauchy stress tensor of classical 


elasticity theory (see also Lazar : (1 20141) ). If one substitutes 
Eq. © into Eq. m, then one obtains 


LCijkl^-kl — & ij 


(13) 


which is a partial differential equation for the elastic 
strain tensor of gradient elasticity of Helmholtz type. 
Eq. © should not be confused with a kind of modi¬ 
fied Hooke’s la w; what i s done in the “basic postulate” 
of GRADELA (lA ifantis. 1992 : Altan and Aifantis . 1992 : 


Gutkin and Aifantisl . 1996 . 1997 ). 

One of the main results of the (LP) article (what (A) 
says as the first major criticism of the (LP)) is that: “Eqs. 
(30) and (31) are not solution for a mode III crack in nonlo¬ 
cal elasticity since they do not satisfy the equilibrium con¬ 
dition (4)”. This is commented by (A) in the framework 
of gradient elasticity and not in the framework of nonlocal 
elasticity causing great confusion to the reader. Actually, 
the comment takes place in the wrong framework. We 
would like to mention again that the non-singular stresses 
of mode III, Eqs. (LP30) and (LP31), as well as of mode I, 
Eqs. (LP47)-(LP49), do not satisfy the equilibrium con¬ 
dition (LP4) of nonlocal elasticity. Thus, the mentioned 
nonlocal stresses cannot be considered as physical solu¬ 
tions in nonlocal elasticity, since they produce forces at 
the corresponding crack tips as shown in Eq. (LP35) for 
mode III and in Eqs. (LP51) and (LP52) for mode I. 

Another important result of the (LP) article (what (A) 
says as the second major criticism of the (LP)) is that the 
non-singular stresses and the corresponding non-singular 
strains of mode III and mode I do not fulfill the com¬ 
patibility conditions in the framework of gradient elastic¬ 
ity and they give rise to non-vanishing dislocation den¬ 
sity contributions at the crack tips. In particular, (A) 
comments the fact that the strains (33) and (34) ex¬ 
hibite d in (LP) article (originally published in lAifantisI 
( 20091) ) are not compatib le by writing that the strains 
(46) in (Lazar and Mauginl (2005) are also not compatible; 
stating moreover erroneou sly that there is a c o ntrad iction 
between (LP) article and Lazar and Mauginl ( 2005 ) con¬ 
cerning this point. It should be here clarified that there is 
no contradiction, since there is a grea t differ ence between 
the two problems. lLazar and Mauginl (|2005r) investigate a 
problem of dislocations in the framework of incompatible 
gradient elasti city theo r y, where th e strains are incompati- 
ble. Whereas, Aifantis ( 20091 2011 ); Isaksson et al.1 ( 20121) 
and llsaksson and Hagglundl (12013 ) investigate a problem 
of cracks and starting from classical compatible elastic dis¬ 
tortions, they obtain incompatible elastic distortions due 
to gradient elasticity. In other words, gradient elasticity 
creates incompatibilities (dislocations) in such a problem 
of cracks as an artifact. 

Furthermore, it is misleading to introduce “microstress” 
and “macrostress” in gradient elasticity as done in (A). 
Concepts of “microstress” and “macrostress” make sense 
in microcontinuum theories like micromorphic elastic¬ 
ity, where microst rain is defined and contained (see, 
e.g., Eringen ( 19991) ). 
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Moreover, we show that the fact that the stresses © and 
© are divergence-free, that is Tijj = 0 and T^k.kj = 0 , is 
an outcome of the theory and it is not restrictive as it is 
claimed by (A) These properties of the s tresses were first 
pointed out bv lLazar and Mauginl ( 2005 !). Let us discuss 
first some properties of the solution of the inhomogeneous 
Helmholtz equation (fT2l) or (LP19) appearing in gradient 
elasticity of Helmholtz type. Using the theory o f partial 
differential equations (see, e.g.. lVladimirov ( 197llU . we ob¬ 


tain the representation of the stress tensor as convolu¬ 
tion of the Green function G of the Helmholtz equation, 

tO 


LG = 6 , with the classical singular stress tensor <t° 




(14) 


which is the (particular) solution of the inhomogeneous 
Helmholtz equation m- 5 is the Dirac delta func¬ 
tion and * deno tes the spatial convolution. As discussed 
in Lazail ( 2014 ). the tensor r,; ? - has the physical mean¬ 


ing of the Cauchy stress tensor of gradient elasticity; an 
interpretation which is in full agreement with t he in- 
terpretation used by Eshel and Rosenfeld ( 1970l Il973l) : 
Gao and Ma ( 2010a bl): Ma and Gao ( 2011 ). By apply¬ 
ing the Helmholtz operat or L to both s ides of Eq. ED, 
the result reads (see, e.g.. iKanwa 3 (SHI)) 


Lth = L(G * a%) = (LG) * = 8 * <r? = a* 


(15) 


which proves that Eq. ED solves Eq. ED- Such solution 
is unique in the class of generalized functions. If we use 
Eq. (ED- the property of the differentiation of a convo- 
lutio n and that the operation of convo lution is commuta¬ 
tive (IVladimirovL 197lP iKanwal . 1983 ). then we find that 
the divergence of the Cauchy stress tensor of gradient elas¬ 
ticity is zero 




13,3 


= (G*cr? i )j = G*(<r9 


i W ) = 0 , (16) 

since of- ■ = 0 (see Eq. (fTTM . Hence, = 0. Next, we 
show that t ijk.kj = 0. In order to differentiate a con- 
volut ion, it suffices to differentiate any one of the fac¬ 
tors ( Kanwall Il983h . Therefore, if the convolution (flTl) 
exists, then the Cauchy stress tensor of gradient elasticity 
is self-equilibr ated. If one u ses a so-called “Bifield” ansatz 
for the stress ( Lazail .2014) 


Tij = cr°- + alj , (17) 

where the tensor <r?„ is the classical stress tensor and the 

L 3 

tensor cjL, which is the gradient part of the str e ss, co r- 
responds to the relative stress tensor (see Lazar ( 2014IH . 


then ajjj = 0 follows from Eq. ED- Moreover, the double 
stress tensor (© (see Eq. (LP11) ) is e quilib rated by the 
relative stress tensor crL (see also Vardoulakis and Sulem 
(1 19951) ) 


®ij — — £ A. 


It is easy to see that the following relation holds 

a ij,j = T ijk,kj = (- = 0 , 


(18) 


(19) 


but t ijk.k 7 ^ 0. Therefore, we have shown that the stresses 
Tij and Tijk are divergence-free, T t jj = 0 and Tijk,kj = 0 , 
and this result arises naturally from the theory. Note 
that the non-singular gradient-enric hed stresses o f dis¬ 
l ocations an d disclinations given by iGutkin and Aifantisl 
( 1999 1 l2000l) are also divergence-free. Additionally, using 
Eq. (ED, the boundary conditions of g radient e lastic ity in 
terms of q£ ; and <j) :j can be found in lArava 1 (l 201 lh and 
iza3(l20ll . 

The r -3 / 2 behavior of the “total” stresses at the 
crack tip in the framework of strain gradient elastic- 


Aravas and Giannakoooulos (20091); Fanniiang et al. 

(2002 

): iGeorgiadis (20031); Gourgiotis and Georgiadis 

(2009 

); Gourgiotis et al. (2010); Shi et al. (|200d); 

Stamoulis and Giannakopoulos ( 

20080) and numerically 

(e.g., Akaranu and Zbib (2006 

); Karlis et alj |2007l); 

Markolcfas et al. (2008); Pananicolonulos et al. (20101)). 


Most of these papers shoul d be known to Aifantis for 
many years (see, e.g., in Aifantisl ( 20091 201ll) L An ap¬ 


propriate treatment of a crack problem in the framework 
of strain gradient elasticity is to solve the homogeneous 
Helmholtz-Navier equation for the displacement vector 
(Eq. (flOll ). since a crack problem is a compatible elasticity 
problem. Moreover, Eq. ED has to be accompanied 
by the corresponding boundary conditions of gradient 
elasticity, the so-called Mindlin boundary conditions (see 
Eqs. (LP10)) which are stemming from variational prin¬ 
ciples, and by the appropriate re gularity conditions . One 
can see clearly this procedure in iGeorgiadia (2003) for a 
mode III crack problem. Then, the “total” stress becomes 
singular with r -3 / 2 singularity at the crack tip, as it was 
predicted in the aforem entioned t h eoretica l and n u merica l 
works. In addition, iGourgiotis and Gcorgiadid (|2009h 
solved the homogeneous Helmholtz-Navier equation ED 
of a compatible elasticity problem for mode I and mode 
II cracks together with the appropriate boundary and 
regularity conditions without using the “Ru-Aifantis 
theorem”. They found that even if the “total” stress 
is hypersingular, the “monopolar” stress tensor which 
is directly connected with the compatible elastic strain 
tensor (and corresponds to the Cauchy stress tensor of 
gradient elasticity in our notation) is bounded in the 
vicinity of the crack tip. Moreover, the strain field is also 
bounded at the crack tip region. 

Later, (A) writes about his students and very personal 
views. In the rest of his comment (A) discusses some un¬ 
published results which are related to forthcoming papers 
and not to the results of (A) criticized in the (LP) article. 
Therefore, there is no need to discuss these items. 


3. Conclusion 

The criticisms of (A) against results in (LP) are re¬ 
futed by means of clear and straightforward arguments 
from the mathematical and physical point of view. The 
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present discussion reveals also that consistency is of ma¬ 
jor significance, since the two comments of (A) against 
the two major results of the (LP) article take place in the 
wrong framework; confusing and mixing nonlocal elastic¬ 
ity with gradient elasticity and compatible elasticity with 
incompatible elasticity. A problem of cracks has a different 
treatment than a problem of dislocations. 

Therefore, we conc l ude again that the non-singular 
stresses given by llsaksson et al. (2012) cannot be con¬ 
sidered as physical crack solutions in nonl ocal elastic¬ 
ity a nd th e non-singular strains g i ven by Aifantis ( 20091. 
2 01 it ) and llsaksson and Hagglundl ( 2013 ) cannot be con¬ 


sidered as physical crack solutions in strain gradient elas¬ 
ticity. In nonlocal elasticity the equilibrium condition 
is not satisfied, and in gradient elasticity the compat¬ 
ibility condition is not fulfilled. In addition, the non ¬ 
singular solutions of Aifantlsl d20111) : llsaksson et al. ( 2012 ): 
llsaksson and Hagglundl (|2013l ) for a mode I crack do not 
satisfy the zero stress boundary condition at the crack 
faces. 
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